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ORIGINAL ARTICLE 

A mesenchymal glioma stem cell profile is related to clinical 
outcome 

A Balbous^'^'l U Cortes^'^'^ K Guilloteau^'^'^ C Villalva^'^'^ S Flamant'^'^ A Gaillard^'^ S Miiin^ M Wager^ N Sorel^'^'^ 
J Guilhot^° A Bennaceur-Griscelli'^'^ ATurhan'^'^ J-C Chomel^'^'^ and L Karayan-Tapon^'^'^ 

Recent studies have demonstrated a relationship between the expression of stem cell-associated genes and relapses in 
glioblastoma (GBM), suggesting a key role for tumor stem cells in this process. Although there is increasing interest in this field, 
glioma stem cells (GSCs) are still poorly characterized, their 'stemness' state and factors maintaining these properties remain largely 
unknown. We performed an expression profiling analysis of pluripotency in gliomaspheres derived from 1 1 patients. Comparative 
analysis between GSCs and HI and H9 human embryonic stem cells as well as H9-derived neural stem cells indicates major 
variations in gene expression of pluripotency factors Nanog and 0CT4, but a stable pattern for 50X2 suggesting its important 
function in maintaining pluripotency in GSCs. Our results also showed that all GSC lines have the capacity to commit to neural 
differentiation and express mesenchymal or endothelial differentiation markers. In addition, hierarchical clustering analysis 
revealed two groups of GSCs reflecting their heterogeneity and identified COLlAl and IFITMl as the most discriminating genes. 
Similar patterns have been observed in tumors from which gliomaspheres have been established. To determine whether this 
heterogeneity could be clinically relevant, the expression of both genes was further analyzed in an independent cohort of 
30 patients with GBM and revealed strong correlation with overall survival. In vitro silencing of COLlAl and IFTMl confirmed the 
effect of these mesenchymal-associated genes on cell invasion and gliomasphere initiation. Our results indicate that C0L1A1 
and IFITMl genes could be considered for use in stratifying patients with GBM into subgroups for risk of recurrence at diagnosis, 
as well as for prognostic and therapeutic evolution. 
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INTRODUCTION 

Glioblastoma (GBM), the highest grade of glioma tumor (WHO 
grade IV), is the most frequent and aggressive primary brain 
tumor in adults. The current standard of care for patients with 
GBM includes tumor resection followed by adjuvant radiation 
therapy and chemotherapy. A significant survival benefit was 
reported for GBM patients treated with temozolomide combined 
with radiation therapy.^ However, despite advances in surgical 
and medical neuro-oncology, the prognosis for GBM patients 
remains poor and median survival is < 1 year. As a result of its 
infiltrating characteristics, complete resection of GBM is not 
possible and recurrence occurs almost invariably at the primary 
location of the tumor. Recent data suggest that tumor relapses 
may be linked to the presence within the tumor of a cell sub- 
population with stem cell characteristics, called 'tumor-initiating 
stem cells' or 'cancer stem cells' (CSCs). Such cells, designated as 
glioma stem cells (GSCs) have been isolated from human glioma 
tissues.^'^ First reports suggested that tumorigenic cells in GBM 
were restricted to the CD133+ population.^'^ However, recent 



studies have revealed that CD133-negative cells isolated 
from GBM could also be tumorigenic^"^ and may in fact 
represent inter-convertible phenotypic states of the same cell 
population giving rise to a heterogeneous GSC population.^ These 
cells are highly resistant to conventional chemotherapeutic drugs, 
including temozolomide^^'" and mediate tumor recurrence 
following radiation therapy.^ ^ Therefore, complete eradication of 
GSCs may be a prerequisite for successful therapeutic strategies. In 
order to arrive at such an approach, in vivo and in vitro models 
that faithfully recapitulate the stem cell compartment of gliomas 
have been developed. Among these models, gliomaspheres, 
also termed neurospheres, are cultured in vitro and enriched in 
GSCs. This in vitro model reflects biological and pathological 
characteristics of the stem cell compartment. Recently, using a 
comparative genomic analysis between GSCs and human neural 
stem cells (NSCs), Sandberg et al.^^ identified, alterations of Wnt 
and Hedgehog signaling pathways as important events conferring 
self-renewal potential and tumorigenic properties to GSCs. 
However, the characterization of molecular events underlying 
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GSCs tumorigenesis remains challenging as GBM tumors 
comprise heterogeneous populations of GSCs with distinct 
molecular alterations associated with contrasting phenotypes 
in tumor initiation and differential response to therapy.^^'^^ 
Human embryonic stem cells (hESCs) represent archetypal stem 
cells that maintain their self-renewal and pluripotency by 
blocking differentiation through a well-coordinated trans- 
criptional network.^ These biological features are shared by 
both normal and CSCs.^^"^^ The ESC-like transcriptional program 
was shown to be activated in poorly differentiated tumors 
including GBMs and appeared to be predictive of metastasis 
and death. ^^'^^ However, unlike ESCs, the stemness abilities of 
different types of CSCs may not be reflected in a single-gene 
expression program and there are likely to be many routes to 
stemness as suggested by Muller et al?^ Stem cell traits normally 
present in human NSCs are shown to be aberrantly co-expressed 
with pluripotent markers of hESCs in gliomas. This general 
deregulated expression profile correlates with increased 
malignancy and probably reflects its cell of origin.^^ To date, 
pluripotency/stemness characteristics, and differentiation 
potential of GSCs remain largely unexplored. A better under- 
standing of these aspects would be crucial for further 
improvements in GBM therapy. To address these questions, we 
took advantage of the unique pluripotential characteristic of 



human embryonic and NSC lines, and compared the gene 
expression profiles of 11 GSC lines isolated from GBM patients 
with those of HI and H9-ESCs, and with HQ-derived NSCs. 

RESULTS 

TaqMan stem cell pluripotency assays in GSCs 
To highlight the pluripotency/stemness status and the differentia- 
tion ability of gliomaspheres, which are considered to be enriched 
in GSCs and represent tumoral stem cells, the expression levels of 
selected genes were compared with those of hESCs (HI and 
H9-ESCS) and human NSCs (H9-HNSCs). The pluripotency and 
stemness status of H9-HNSCs, which serves as a reference in our 
strategy, was first compared with the well- characterized HI and 
H9-ESCS. We observed a lower expression of Nanog and P0U5F1 
in H9-HNSCS, but a moderately elevated expression of 50X2 
(Supplementary Figure 1A). As expected, NSC and neural 
morphogenesis genes {Nestin, NeuroDl and PAX6) as well as 
GCMl (glial cells missing gene) having a neurogenic role, were 
expressed at a higher level in H9-HNSCs. On the other hand, glial 
fibrillary acidic protein (GFAP) (highly specific for cells of astroglial 
lineage), oligodendrocyte and neuron development transcription 
factors {0LIG2, HLXB9 and ISLI) were found to be equally 
expressed or downregulated (Supplementary Figure IB). In 
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Figure 1. TaqMan low-density array analysis, (a) mRNA expression of stemness and pluripotency-related genes in GSCs as compared with the 
expression of H9-HNSC. (b) mRNA expression of stemness and pluripotency-related genes in GSCs as compared with the expression of HI, 
H9-ESC lines. (c,d) Expression analysis of neural (GFAP, NEUR0D1 and 0LIG2) endothelial {CD34 and PECAM1) and trophoblastic {GCMl) 
differentiation markers in GSCs as compared with H9-HNSCs. The chart shows median values with range from 1 1 GSC lines determined by 
2-AAct reiai-jvg quantification method. Relative expression of target genes determined by the using GAPDH as endogenous control and 
normalized to H9-HNSC or HI, H9-ESCs calibrator. 
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addition, our results confirmed that H9-HNSCs are multipotent 
because the genes involved in mesoderm (for example, C0L1A1, 
RUNX2, C0L2A1 and CDH5) or in endoderm development (FA/7, 
MM/A 7, LAMBl and LAMCl) have a relative high constitutive 
expression (Supplementary Figure 1C). 

The pluripotency and stemness status of GSCs was next 
compared with that of H9-HNSCs cultured under the same 
conditions. Nanog and P0U5F1 were found to be equally 
expressed in H9-HNSCs and GSC lines, whereas the expression 
of 50X2 appeared to be lower in GSCs (Figure 1a). Interestingly, 
the expression of 50X2 in all GSCs was quite similar to H1 
and H9-ESCS (<1 log; Figure lb). Nestin and GCM7 were 
also found to be equally expressed in H9-HNSCs and GSCs 
lines while NeuroDl and PAX6 expression were lower in GSCs. 
In contrast, genes involved in astrocyte (GFAP), neuron 
(HLXB9) and oligodendrocyte (0UG2) differentiation were found 
to have a higher expression in GSCs indicating that GSCs are 
more committed to neural differentiation than H9-HNSCs 
(Figure 1c). However, GSC population showed multipotent 
potential because genes involved in endothelial differentiation 
such as CD34 or PECAMl had expression patterns similar to 
H9-HNSCS and the endoderm-related transcription factors 
GATA4 were found to be overexpressed in GSCs compared with 
NSCs (Figure Id). 

A hierarchical clustering analysis was performed using Ward's 
linkage method and Euclidean distance. The unsupervised 



method segregated the 1 1 GSCs samples from H9-HNSCs, 
HI- and H9-ESCs, and clearly separated GSCs into two distinct 
groups (Figure 2a). Group GSCX-1 included GSCs 1, 2, 4, 6, 8, 10 
and group GSCX-2, GSCs 3, 5, 7, 9, 11. To identify genes 
differentially expressed between the two groups, individual mRNA 
expression levels were compared (Supplementary Table 3). Seven 
genes [COLlAh IFITMh SFRP2, FNl, 557, F0XA2 and GAL) presented 
a significant higher expression in group GSCX-2 as compared with 
group GSCX-1 using a f-test with Benjamini-Hochberg correction 
(Supplementary Figure 2). It is noticeable that C0L1A1 and IFITMI 
expression harbored the most significant difference between the 
two groups of patients (Figure 2b). The expression level of these 
two genes was also performed in bulk tumors from which GSCs 
were established. Interestingly, the expression of both genes was 
shown to be higher in bulk tumors from group TX-2 than in bulk 
tumors of group TX-1 (Figure 2c). However, no significant 
difference in overall survival was observed between both groups 
of patients (Table 1). 



Effects of C0L1A1 and IFITMI inhibition on cell invasion, 
proliferation and neurosphere initiation 

In order to determine the biological significance of differential 
expression levels of COLlAl and IFITMI in GSCs, we inhibited both 
gene expression and studied their subsequent effect on invasion, 
proliferation and neurosphere initiation. Lentiviral vectors 
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Figure 2. Hierarchical clustering analysis of GSCs. (a) The unsupervised clustering using Ward's linkage method and Euclidean distance 
individualized H9 HNSCs, H1-H9-ESCs and two groups of GSCs (GSCX-1 and GSCX-2). Note that culture conditions of H9 HNSC do not 
influence the expression of pluripotency gene samples H9 HNSC-1 and H9 HNSC-2 samples, (b) Box plot representation of C0L1A1 and /F/rM7 
expression levels between GSC groups GSCX-1 and GSCX-2. The charts show the logio expression of relative quantification (RQ) values 
normalized to the expression of H9 HNSC cell line, (c) Box plot representation of COLIAI and IFITMI expression levels between bulk tumor 
groups TX-1 and TX-2. The top edge of the boxes represents the 75th percentile, the bottom edge, the 25th percentile and x the mean. The 
range is shown as vertical edge. 
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Table 1. 


Characteristics of GBM patients from whom neurospheres were derived 
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Abbreviations: GBM, glioblastoma; GSC, glioma stem cell; RX, radiotherapy 60 Gy; TMZ, temozolomide. 





Figure 3. Knockdown efficiency of C0L1 Al and IFITMI and effect on cell invasion, (a) C0L1 Al and IFITMI nnRNA expression was nnonitored by 
quantitative reverse transcriptase-PCR in GSC-3 and GSC-9 {n = 2). Significant inhibition of 90% of nnRNA expression was observed. GADPH 
was used as internal control. *P<0.05. Each value is expressed as nnean ± s.e.nn. (b) Western blot analysis showed dinninution of IFITMI and 
C0L1A1 proteins expression, p-actin was used as loading control, (c) Densitonnetry of protein level. The bar graph shows relative IFITMI and 
C0L1A1 protein expression levels expressed as percentage of IFITMI and C0L1A1 expression in control cells {n = 3). We observed >80% 
inhibition for both genes in GSC-3 and GSC-9 cells, as connpared with control cells, (d) Representative fields of Boyden cell invasion assay. After 
3 days of invasion, cells that nnigrated across the mennbrane were fixed and stained with hennatoxylin/eosin and were counted in one hundred 
different fields. The invaded cells were photographed under the nnicroscope at x 100 nnagnification. (e) A significant reduction of invasion is 
observed in both GSC-3 and GSC-9 cells lines with sh-IFITMl as well as for sh-COLIAI. The results represent a percentage of invaded cells. 
Data are presented as nnean ± s.e.nn. fronn three independent experinnents. ^P<0.05, ^^P <0.01. 



encoding short hairpin RNA (shRNA) directed against COLlAl and 
IFITMI were used to transduce GSC-3 and GSC-9, which displayed 
a high basal expression of the two genes. Inhibition of COLlAl 
and IFITMI expression were nnonitored by quantitative reverse 



transcriptase-PCR and western blot. For the two genes, a decrease 
of >90% in nnRNA expression was observed in GSC-3 and GSC-9 
as compared with control cells (P<0.05; Figure 3a). Sinnilarly, 
western blot analysis revealed an 80% reduction for C0L1A1 and 
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IFITMI proteins levels in GSC-3 and GSC-9 cells as compared with 
control cells (Figures 3b and c). 

The effect of shRNA-C0/./>A7 and shRH A- IFITMI on cell invasion 
was first studied using Boyden cell invasion assay. After 3 days, 
we observed a significant reduction in the invasive potency of 
sh-IFITM1 GSC-3 (72%) and sh-IFITMl GSC-9 (67%) cells as 
compared with control cells (P<0.05). An equivalent reduction 
of invasion was also shown for sh-COLIAI GSC-3 (72%) and 
sh-C0L1A1 GSC-9 (70%) cells (Figures 3d and e). The effect of 
inhibiting COLIA and IFITMI expression on cell proliferation was 
then analyzed using a MTS assay. GSC-3 and GSC-9 cells have a 
6-day and 1-day doubling time, respectively. Inhibition of IFITMI 
or COLlAl expression had no effect on doubling time for either 
GSC lines (5.7 days for sh-IFITMl GSC-3, 2 days for sh-IFITMl GSC-9 
and 4.6 days for sh-C0L1A1 GSC-3). Nevertheless, we observed a 
slight increase of doubling time for sh-C0L1A1 GSC-9 (2 days; 
Figure 4a). Finally, the effects of inhibiting COLlAl and IFITMI 
expression on neurosphere-forming capacity were evaluated. The 
frequency of neurosphere-initiating cells corresponding to the 
number of cells required to generate at least one tumor sphere 
per well was 1/25 for GSC-3 and 1/15 for GSC-9. Inhibition of 
IFITMI expression only slightly modified the neurosphere-initiat- 
ing capacity of GSC-3 (1/27) and GSC-9 (1/21). However, inhibition 
of COLlAl resulted in a decrease of neurosphere-initiating 
cell frequency, as the number of cells required to form neuro- 
spheres was 1/128 for sh-COLIAI GSC-3 and 1/35 for sh-COLIAI 
GSC-9 (Figure 4b). 

Altogether, these observations indicate that C0L1 Al and IFITMI 
knockdown markedly decrease the invasive capacity of GSCs 
from primary human GBMs. Furthermore, suppressing COLlAl 
expression strongly affects the neurosphere-initiating capacity 
of GSCs. 



Expression of COLlAl and IFITMI in tumors and relation to survival 
The absence of significant difference in overall survival 
betweenTX-1 and TX-2 groups may be due to the small number 
of cases. To test this hypothesis and to determine whether the 
expression level of COLlAl and IFITMI in tumors could be of 
clinical significance, we analyzed the mRNA expression of COLlAl 
and IFITMI in primary tumor samples of an independent cohort of 
30 GBM patients (Table 2). A hierarchical clustering separates 
patient samples into two statistically distinct subgroups (called A 
and B). Patients from subgroup A (n = 13) showed a significant 
reduction in the expression of COLlAl and IFITMI as compared 
with patients from subgroup B (P = 0.037 and P = 0.022, 
respectively; Figure 5a and Supplementary Table 4). Kaplan- 
Meier survival analysis demonstrated that patients from subgroup 
A had a median overall survival of 10.5 months (95% confidence 
interval: 8.2-14.0) and patients from subgroup B had an overall 
survival of 17.5 months (95% confidence interval: 15.5-24.7; 
P<10"'^; Figure 5b). 



DISCUSSION 

Recent data suggest that GBMs belong to a tumor class in which 
the CSC hypothesis may be valid with the presence of a small 
population of cells displaying self-renewal and differentiation 
abilities. These progenitors have been detected by the use of 
In vitro and in vivo transplantation assays^"^ with a hierarchy 
similar to that found in hematopoietic stem cells.^^ Several 
malignant features of GBMs such as recurrence, radio- and chemo- 
resistance may be due to the presence of GSC population, which 
could represent a therapeutic target. The purpose of the present 
work was to investigate the pluripotency, stemness characteristics 



a GSC-3 b GSC-3 




Figure 4. Effects of C0L1A1 and IFITMI inhibition on cell proliferation and on neurosphere-initiating capacity, (a) Cell proliferation was 
measured by using a MTS assay. Each experiment was carried out in replicates of six and repeated twice and expressed as mean ± s.d. 
(b) Inhibition of COLlAl resulted in a significant (*P<0.05) decrease of neurosphere-initiating cell (NS-IC) frequency whereas inhibition of IFITMI 
expression only slightly modified the neurosphere-initiating capacity. Final cell dilutions yielding 37% of negative wells correspond to the 
dilution at which there is one NS-ICs per well. Each experiment was performed in duplicate and repeated twice and expressed as mean ± s.d. 
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Table 2. Characteristics of 30 GBM patients from an independent cohort study 
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and the differentiation potential of GSCs using a Taqnnan Low 
Density Array (TLDA)-based connparative strategy. Our results have 
highlighted a sizable decrease in Nanog and P0U5F1 nnRNA 
expression in hNSCs as well as in GSCs connpared with pluripotent 
hESCs. These data are in line with those of a previous study from 
Choi et al}^ showing, through cell reprogrannming experiments, 
that these two genes had a key role in controlling the pluripotency 
and the differentiation status of hNSCs. Among genes involved in 
the maintenance of pluripotency or expressed in undifferentiated 
cells, only 50X2 showed a similar expression in GSCs and H9-ESCs. 
This observation suggests that 50X2 largely contributes to the 
maintenance of a pluripotent and undifferentiated state in GSCs. 
50X2 is known to be involved in NSC homeostasis and its 
expression is essential for the maintenance of induced pluripotent 
stem cells as well as long-term in vitro neurosphere cultures.^^'^° 
On the basis of TLDA results, the differentiation potential of 
GSCs was consistent with their cellular origin. Genes involved in 
astrocytic [GFAP), neuronal [NFS and NEUR0D1) and oligodendro- 
cytic (0LIG2) differentiation were clearly overexpressed as 
compared with hNSC. Similarly, the GCM7 gene, which is thought 
to act as a binary switch between neuronal and glial cell 
determination, displayed higher expression in GSCs. In the 
Drosophila nervous system, the ortholog of GCM7 is transiently 
expressed in glial precursors to switch their fate from the 
neuronal default state to glia.^^ Interestingly, the genes involved 
in bone {RUNX2) and endothelial {CDH5, CD34 and PECAM) 
differentiations are also overexpressed in hNSCs and GSCs 



as compared with HI and H9-ES cells. Plasticity is a 
characteristic of NSCs and differentiation of NSCs and GSCs 
into unanticipated cell types such as endothelial cells has already 
been shown.^^"^^ 

Hierarchical clustering based on TLDA analysis clearly identified 
two GSC groups. Seven genes appeared significantly over- 
expressed in group GSCX-2. Among them, COLlAl, FNl and 
LAMB! known to be mesenchymal markers are involved in the 
extracellular matrix. The role of extracellular matrix components 
such as fibronectin (FNl) in glioma invasion is now well known. In 
this regard, microarray studies have revealed an overexpression of 
FA/7 in human gliomas as compared with normal brain tissue.^^ A 
recent study performed by Kim et al.,^^ which aimed at identifying 
predictive survival models, classified GBM into three distinctive 
subgroups according to the differential expression of 
mesenchymal-related genes. Interestingly, the mesenchymal- 
enriched subgroup of GBM was found to significantly 
overexpress COLlAl and FNl genes, as observed in our study in 
group GSCX-2. Laminin subunit (31 (LAMB1) was also found 
overexpressed in human gliomas. Among genes differentially 
expressed in the two GSC groups, C0L1A1 and IFITMI were found 
to be the most discriminating. COLlAl is considered as a 
mesenchymal osteoblast marker^^ and was also defined as a 
glioma endothelial marker selectively expressed in the 
microvasculature.^^ In vitro silencing of C0L1A1 in two different 
GSC lines revealed its prominent role in promoting invasion and 
neurosphere-initiating capacity. IFITMI is a transmembrane 
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Figure 5. C0L1A1 and IFITMI expression levels in an independent cohort of 30 patients' tumors, (a) Box plot representation of the mRNA 
expression levels. The charts show the logio expression of relative quantification (RQ) values normalized to the expression of H9 NSC cell line. 
The top edge of the boxes represents the 75th percentile, the bottom edge, the 25th percentile and x the mean. The range is shown as 
vertical edge, (b) Kaplan-Meier analysis for overall survival of an independent cohort of 30 GBM patients. Survival distribution of patients with 
low (group A; n = 13) versus high (group B; n = 17) expression of C0UA1 and IFITMI (P<10"'^). 



protein whose expression is induced by interferon y. It is 
expressed on the cell surface of primordial germ cells in a 
developmentally regulated manner and is supposed to 
modulate cell adhesion and differentiation.^^ This protein is also 
involved in the invasion process and progression of gastric 
cancer cells,'^° colorectal tumors"^^ and head and neck squamous 
cell carcinomas'^^ through the activation of matrix 
metalloproteases -12 and -13. Recently, gene invalidation 
studies using shRNA unveiled the role of IFITMI in proliferation, 
migration and invasion of GBM cell line U87MG.^^ Our 
results obtained In vitro using two IFITMI knockdown GSC 
lines are corroborate previous observations because IFITMI 
silencing led to a significant reduction of invasive capacity. 
Taken together, our results highlight two distinct types of 
glioma-initiating stem cells according to the expression 
of mesenchymal markers. 

The expression level of C0L1A1 and IFITMI in bulk tumors from 
which the GSCs were established probably reflects their expres- 
sion level in GSCs. The lack of difference in overall survival 
between the two groups is probably due to the small number of 
cases (n = 11). However, analysis of an independent cohort of 30 
primary GBM patients clearly identified two groups of patients 
according to the expression levels of COLlAl and IFITMI in their 
tumors with a statistically significant difference in overall survival, 
which increased from 10 to 17.5 months. Importantly, the groups 



of patients were not different in terms of age, ECOG score or 
treatment regimens consisting in concomitant temozolomide and 
RT after surgery. Further, no alteration in factors known to 
influence patients prognosis {IDHI and IDH2 mutations or MGMT 
methylation status)^^'^°''^^''^^ was observed in either group. 
In addition, class comparison between very long-term survivors 
(>24 months) and very short-term survivors (<12 months) in the 
Rembrandt database revealed a significant clinical value for 
COLlAl (P = 0.04), in agreement with our results (https:// 
caintegrator.nci.nih.gov/rembrandt/). 

A recently established notion postulates that neoplastic 
transformation in the central nervous system converts neural 
cells into cell types manifesting a mesenchymal phenotype 
state associated with uncontrolled ability to invade and stimulate 
angiogenesis.'^^''^'' Thus, a mesenchymal phenotype would be 
the hallmark of tumor aggressiveness in human malignant 
glioma in contrast to proneural phenotype associated with 
longer survival.^^'"^^ Consistent with this concept, recent reports 
have suggested an epithelial-to-mesenchymal-like transition 
phenomenon in GBM, similar to epithelial tumors but lacking 
cadherin reorganization."^^ 

In a recent in vitro study, Mao et aC^ have characterized 
mesenchymal GSCs using a 5796-gene signature and reported 
an increase in glycolytic pathways in these cells. These authors 
also showed promising therapeutic opportunities through 
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inhibition of aldehyde dehydrogenase 1 A3 in these 
mesenchymal-enriched GSCs. In our study, overexpression of 
IFITMI and COLlAl appears predictive of a longer survival. This 
result may seem paradoxical as it is known that GBMs with 
mesenchymal phenotype are aggressive tumors. However, recent 
data have shown a beneficial effect with better prognostic in 
patients with the mesenchymal subtype,^^ supporting the 
difference of survival observed in our study for patients with 
high levels of COLIAl and IFITMI expression. In conclusion, the 
expression of these two genes could represent useful molecular 
markers, not only to characterize GSCs, but also to stratify 
GBM patients at diagnosis since both genes appears to correlate 
with a longer survival. 

MATERIALS AND METHODS 

Patients 

Gliomaspheres were derived from 1 1 adult patients with high-grade 
gliomas (Table 1). Tissues were obtained during surgery (Department of 
Neurological Surgery, Poitiers University Hospital, France) after informed 
consent of all patients and in accordance with the Declaration of Helsinki. 
Patients were free from any past medical history, especially with regard to 
cancer. Tumor diagnosis and grading were established according to the 
WHO criteria^^ and were systematically confirmed by two expert 
neuropathologists. With the approval of the ethics committee, validation 
of COUAl and IFITMI expression was performed in a panel of primary 
tumors from an independent cohort of 30 patients (Table 2). 

Determination of IDHI, IDH2 and MGMT promoter status in GBM 
patients 

The genomic region spanning wild-type R132 of IDHI and R172 of IDH2 
was analyzed by direct sequencing as previously described.^^ Methylation 
status of the MGMT promoter was determined by bisulfite modification 
and subsequent nested methylation specific-polymerase chain reaction.^^ 

In vitro and In vivo GSCs characterization 

Tumor samples were obtained within 30min after surgical resection from 
11 adult GBM patients (GSC-1-GSC-1 1) and tumor sphere cultures were 
performed and characterized as we described previously.^°'^^ Isolation and 
characterization of GSCs were based upon criteria described by Singh 
et al} and Galli et ai^ All GSC lines were assessed for in vitro clonogenicity 
by limiting dilution assays. The number of cells required to form one 
neurosphere, which reflects the frequency of CSCs in the entire population, 
was determined as previously described^° (Supplementary Table 1). 

hESC culture and H9-derived human NSC culture 
After authorization from the French Biomedicine Agency (permission to 
conduct research on human embryonic stem cells, SANB0721164S and 
SANB0721165S, 9 February 2007), hESC lines HI and H9 (WiCell Research 
Institute, Madison, Wl, USA) were cultured on mitomycin-treated mouse 
embryonic fibroblast cells in KnockOut Dulbecco's modified Eagle medium 
supplemented with 20% KnockOut serum replacement (Invitrogen, 
Carlsbad, CA, USA), and 5 ng/ml basic fibroblast growth factor. GIBCO 
Human Neural Stem Cells (H9-HNSCs) are derived from NIH-approved H9 
(WA09) hESCs (Life Technologies, Grand Island, NY, USA). Cells were 
cultured following the manufacturer's instructions (H9-HNSC-1) as well as 
in neurobasal glioma stem cells medium (H9-HNSC-2). Culture conditions 
of H9-HNSCS did not influence the expression of pluripotency genes 
(Figure 2a). 

Analysis of gliomaspheres mRNA by human stem cell pluripotency 
low-density array 

TaqMan Human Stem Cell Pluripotency Array (Applied Biosystems, Foster 
City, CA, USA) consisted in 4 identical 96 gene sets (90 target genes and 6 
endogenous controls). The list of target genes was based on a previous 
work from the International Stem Cell Initiative^^ and is detailed in 
Supplementary Table 2. These genes were selected with the following 
criteria: (i) expression in undifferentiated cells (embryonic carcinoma cells, 
primordial germ cells or hESCs), (ii) maintenance of pluripotency, 
(iii) expression level correlated with stem cell state (stemness) and 



(iv) determination of the differentiation capability into all three embryonic 
germ layers (ectoderm, endoderm and mesoderm). Briefly, 2 jig of total 
RNA obtained from cells at passage 6 were reverse transcribed using the 
High Capacity RNA-to-cDNA Kit according to the manufacturer's instruc- 
tions (Applied Biosystems). The genetic stability of early passage GSCs in 
serum-free conditions has previously been demonstrated.^^'^^ Each line of 
the microfluidic card was loaded with 100|il of a mixture containing the 
complementary DNA template and 50|il of 2X TaqMan Universal PCR 
Master Mix. Quantitative PCR was then carried out with the ABI PRISM 
7900HT Sequence Detection System (Applied Biosystems). Each 
experiment was conducted in triplicate and target genes that showed 
no amplification, abnormal amplification curves or bad replicates were 
excluded (Supplementary Table 2). Amplification data were analyzed using 
the StatMiner 3.0 software (Integromics, Madrid, Spain). Relative quanti- 
fication of target gene expression was determined by the 2~^^*^^ method 
using glyceraldehyde 3-phosphate dehydrogenase (GAPDH, most stable 
reference gene) as an endogenous control and normalized to the 
expression in H1-H9-ESCs or H9-HNSCs. 



Vector design and preparation of lentiviral supernatant 
shRNA constucts against human COLIAl and IFITMI were cloned 
downstream to human HI promoter into a modified pTRIP lentiviral 
backbone carrying the selectable cassette GFP-2A-puroR under the control 
of the human phosphoglycerate kinase promoter. Small interfering RNA 
sequences and details about the constructs are available on request. 
Lentiviral particles were produced using standard methods. Briefly, viral 
plasmid was introduced together with packaging plasmids into 293T cells 
using calcium phosphate transfection procedure. 293T cells are grown in 
Dulbecco's modified Eagle medium supplemented with 10% fetal calf 
serum and lOOIU/ml of penicillin-streptomycin (Life Technologies). 
Medium was replaced after 24 h of transfection, and viral particles were 
collected after 48 h by ultracentrifugation at 25 000r.p.m. for 2h. Viral 
particles were re-suspended in serum-free medium, aliquoted and stored 
at -80°C until use. GSC-3 or GSC-9 cell lines (1x10^ cells) were 
incubated with the indicated lentiviral particles at a multiplicity of infection 
of 1 for 48 h. Control cells were infected with an empty vector. Two days 
after transduction, cells were selected with 1 jiivi of puromycin (Sigma, 
Lyon, France). 



Cell proliferation 

Doubling times of gliomaspheres as well as the inhibitory effect of shRNA- 
C0L1A1 and shRNA-IFITMl on their proliferation were assessed by CellTiter 
96 Aqueous Non-Radioactive Cell Proliferation Assay (Promega, Lyon, 
France). Cells were plated in a 96-well plate at a density of 5 x 10^ cells per 
well in 100|il medium and quantification of viable cells was performed at 
492 nm with a micro-plate reader (Dynex Technologies, Chantilly, France). 



Cell invasion assay 

Cell invasion assay was performed using Boyden chamber with 8|im pore 
size (Becton Dickinson Biosciences, Le pont de Claix, France), coated with 
matrigel (Becton Dickinson Biosciences). Control cells, shRNA-silenced cells 
were seeded in the upper chamber of each well and NBE medium was 
added to the lower chamber. Cells were allowed to invade for 3 days and 
the upper side of the membrane containing non-invasive cells was cleaned 
with a cotton swab. Invading cells were fixed with cold methanol followed 
by hematoxylin/eosin staining. The membrane was cut-off and the number 
of migrating cells was counted under a microscope (Nikon Eclipse 801, 
TribVN ICS framework, Chatillon, France). 



Limiting dilutions and neurosphere-initiating cell assay 
To determine the effect of shRNA-COLIAI and shRNA-IFITMl on the 
frequency of 'neurosphere-initiating cell', we performed limiting dilution 
assays using single-cell populations immediately after dissociation of 
neurospheres. Final cell dilutions ranged from 1 cell per well to 80 cells per 
well. At day 21, the fraction of wells free of neurosphere for each cell 
plating density was determined. These results were plotted against the 
number of cells plated per well. The number of cells required to form one 
neurosphere, which reflects the frequency of CSCs in the entire population, 
was determined as described previously.^° 
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Western blot analysis 

Cells were lysed in Laemmli buffer and protein concentration was 
determined using a Bradford assay (Bio-Rad, Marnes-La-Coquette, France). 
Equal amounts of protein samples were separated by sodium dodecyl 
sulfate-polyacrylamide gel electrophoresis and transferred onto a nitro- 
cellulose membrane (Amersham Biosciences, Velizy-Villacoublay, France). 
Membranes were blocked with 5% non-fat milk in phosphate-buffered 
saline 0.1 % Tween for 2 h and incubated overnight at 4 °C with the primary 
antibody. Antibodies used in this study were goat anti-IFITMl antibody 
(R&D Systems Europe, Lille, France), rabbit anti-COLIAI antibody (Abeam, 
Paris, France) and mouse p-actin antibody (Abeam). After incubation with 
the appropriate secondary antibodies (Santa Cruz, Dallas, TX, USA), blots 
were revealed by chemiluminescence (ECL Plus reagent, Amersham 
Biosciences). Band intensity was quantified with the use of ImageJ 
software (Bethesda, MD, USA). 



Determination of C0L1A1 and IFITMI mRNA expression by reverse 
transcriptase-quantitative PGR 

Total RNA was extracted from tumor specimens using Qiagen RNeasy Mini 
Kit according to the manufacturer's instructions (Qiagen, Courtaboeuf, 
France). In all, 3|ig of total RNA was transcribed into complementary DNA 
using Superscript II (Life Technologies). Expression of COLlAl, IFITMI, and 
GAPDH was determined using TaqMan Gene Expression Assays 
Hs00164004_m1, Hs00705137_s1, Hs Hs99999905_m1 (Applied Biosys- 
tems) respectively. Relative quantification was established by the 2"^^*^^ 
method with GAPDH as an endogenous control and normalized to the 
expression of the H9-HNSCs. Analysis of the results was done using the 
StatMiner software. 



Statistical analysis 

Clustering gene expression patterns were determined using hierarchical 
algorithms of StatMiner software. In this study, Euclidean distance and 
Ward's linkage method were applied although clustering results were 
identical whatever the linkage method and distance measure used. In 
order to identify target genes with significant differences in their 
expression, a parametric paired f-test was used adjusting the P-value with 
the Benjamini-Hochberg false discovery rate-controlling procedure to 
obtain reliable significance levels with the high number of tests performed 
(StatMiner). 

An independent cohort of patients (n = 30) was analyzed using SAS v.9 
software (SAS, Gary, NC, USA). Fisher's exact test was used for testing the 
association of groups from cluster analysis with those given by patient's 
characteristics (MGMT). Correlations between gene expression data were 
analyzed using Spearman's rank test. Overall survival estimates with 95% 
confidence intervals were provided by the Kaplan-Meier method and 
compared with the log-rank test. 
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